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Abstract 

We prove localization and probabilistic bounds on the minimum level spacing 
for a random block Anderson model without monotonicity. Using a sequence of 
narrowing energy windows and associated Schur complements, we obtain detailed 
probabilistic information about the microscopic structure of energy levels of the 
Hamiltonian, as well as the support and decay of eigenfunctions. 
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1 Introduction 


1.1 Background 

We present a new method for proving separation of eigenvalnes for disordered qnan- 
tnm mechanical Hamiltonians on the lattice, in the strong disorder or weak hopping 
regime. The traditional approach to this question is based on Minami’s method m 
When available, it provides a bound on the probability of two or more eigenvalues in 
an interval, and thus gives probabilistic control on the separation of eigenvalues. Vari¬ 
ants of the Minami estimate have been an important component of proofs of Poisson 
statistics for eigenvalues in an interval [271[2ni[221[ini[3l[Hl[I7l[I3]. They are important 
for understanding the behavior of the ac conductivity [23]. They also lead to mini¬ 
mum level-spacing estimates and simplicity of the spectrum, provided the probability 
distribution of the potential is Holder continuous with index a > | [21 Ej. 

Obtaining Minami-type estimates has been problematic for so-called ‘non-monotone’ 
models, in which the variation of the single-site potential with the random variables may 
change in sign from point to point in space. While there are a number of results on 
localization for non-monotone models [25] [3ll [H El ES] [IS] [H El [T] [T2|, there are 
only a couple of multi-level results that we are aware of. Simplicity of the spectrum 
in inhnite volume is established in [28] . In [30], a Minami estimate is proven for some 
non-monotone models that admit a transformation to a monotone situation. In m an 
iV-level Minami estimate is proven for a class of random block models, in which & k x k 
random block Hamiltonian sits at each site of the lattice, and the blocks are coupled 
through a deterministic hopping matrix. The randomness must be sufficiently rich to 
ensure the success of the method. Of particular interest is the case of 2 x 2 blocks of 

the form ( ); Minami estimates were obtained in [11] in the situation in which 

\ V —uj ' ’ 

both u and v are /^-regular random variables. Still, this leaves out the case where v is 
deterministic, which may be called a fully diluted model (meaning that there is only 
one degree of freedom in the randomness for each block). 

In this work we consider a fully diluted random block Hamiltonian with 2x2 block 
u 1 \ 

This is a model originally proposed by Spencer; localization and Holder 


1 —u 

continuity of the integrated density of states are proven in [T2|. We prove a probabilistic 
statement on minimum eigenvalue spacing. Specifically, we give a bound showing that 
the probability that the minimum eigenvalue separation is less than 6 tends to zero with 
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S. (See Theorem 3.9 below). Our interest in this problem stems from recent work 
in which a minimum eigenvalue separation condition (called limited level attraction) 
was used as an assumption in a proof of many-body localization for a one-dimensional 
spin chain. Although the separation condition remains unproven in the many-body 
context, one may perhaps gain some insight into the problem by replacing the tensor 
product of 2 X 2 spin spaces with a direct sum, and this leads to the random block 
Hamiltonian considered in this paper. The dilution of randomness is not as extreme as 
in the many-body case, but it remains an important issue nevertheless. 

The method should be of interest in its own right. We use a sequence of Schur 
complements to reduce the Hilbert space to subspaces of smaller and smaller dimension. 
At the same time, we shrink the width of the energy interval under consideration, and 
eventually only one eigenvalue remains. The subspaces are structured in physical space 
as a set of extended blocks whose (renormalized) spectrum is in a small neighborhood 
of some energy E. The Schur complement formula produces effective Hamiltonians for 
the blocks exhibiting exponential decay in the distance between blocks. The blocks 
themselves become more and more dilute as the procedure proceeds, so the distance 
between blocks grows, and hence inter-block matrix elements tend rapidly to zero. The 
key probabilistic input involves estimating the probability that a block has spectrum 
close to E, and bounding the probability that a block has a very small level spacing. 
The determinant and the discriminant of the Schur complement matrix for a block is a 
polynomial in the random variables (the degree depends on the size of the block). This 
ensures a degree of transversality of eigenvalues of subsystems (and their differences) - 
although not uniformly in the size of the subsystem. We obtain a rather weak bound 
on the probability of spectrum near E (replacing the usual Wegner estimate [M]) and 
the probability of near-degeneracy (replacing the usual Minami estimate). Although 
the bound degenerates as the size of the block gets large, large blocks are improbable 
to begin with, and the trade-off determines the form of our results on density of states 
and level separation. There is some similarity with the method of [5], [26], in which 
localization is proven for a Hamiltonian depending analytically on the random variables. 
In those works, the Wegner estimate is replaced with a transversality condition arising 
from polynomial approximation of determinants (using a version of Cartan Lemma 
proven in 0). 

As in 121 . we give a direct construction of eigenfunctions and eigenvalues in terms of 
convergent graphical expansions that depend on the random variables only within each 
graph. However, we gain some precision by working near a hxed energy (as advocated 
by 0)- We have also been influenced by ideas in |3l |2| (the Feshbach map) and of 
course by 


1.2 Model 

We consider the following non-monotone Anderson model. A particle hops in A, a 
rectangle in Z'^, with each position a: G A having two states, -|- and -. The single- 
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hx — 


( 1 . 1 ) 


position Hamiltonian is 

Ux 1 
1 'Ux 

where u = { ux } x&k is a collection of bounded iid random variables with a bounded 
density. For simplicity, let us assume each Ux has a uniform distribution on [-1,1]. The 
Hilbert space is where n = |A| is the number of elements of A. Let Hq = ©xeA^x- 
The hopping matrix is 

1/ = 7 J (g) J 2 , (1.2) 

with I 2 denoting the 2 x 2 identity matrix for the internal space, 7 is a small parameter, 
and 

1 , if lx- 2/1 = 1 ; 

0 , otherwise. 


Jxy 


(1,3) 


Here |x| = Finally, we dehne the Hamiltonian 

H = Ho + H. 


(1.4) 


One could of course consider a much wider class of models. However, this choice will 
serve to illustrate the method in a simple context exhibiting both non-monotonicity 


and dilution. That is, the “bare” energies Ux,—Ux in (1.1) move in opposite directions 


as Ux varies, and there is only one random variable Ux for each 2x2 block of Hq. 


1.3 Main Results 

Our goal is hrst to prove localization by controlling the density of states and by showing 
stretched exponential decay of the eigenfunction correlator \^a{x)ipa{y)\- Second, 

we prove probabilistic bounds on level spacing by constructing individual eigenstates 
and controlling their separation in energy from the others. 


We prove several facts about the random block Hamiltonian (1.4) 


Theorem 1.1. (Theorem 3.1). There exists a constant b > 0 such that for'y sufficiently 


small the following is true. For any 0 < 5 < 7 ^^^, any rectangle A, and any interval 
I = [E — 5/2,E + d/2], let N{I) denote the number of eigenvalues of El in I. Then 


EN(I) < 4|A|exp(-6|log7|^''7ogi5|‘''“). 


(1.5) 


Although this bound goes to zero faster than any power of l/|log5|, it does not 
imply Holder continuity of the integrated density of states. Thus it is weaker than the 
result of [T^ . 


Theorem 1.2. (Theorem 3.1). Let x= g- There exists a constant q> 0 such that for 


7 sufficiently small, the eigenfunction correlator satisfies 


( 1 . 6 ) 


for any rectangle A. 
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A similar bound is proven in [T2] , but with exponential decay. If we restrict the sum 


in ( 1 . 6 ) to eigenvalues in /, we prove a similar bound but with a prefactor tending to 


zero like the right-hand side of (1.5) - Corollary 3.8 We prove these results partly as a 


warm-up to the level-spacing problem - but it is also helpful to understand how strong a 
result is possible with our method and to make comparisons with the fractional-moment 
results. 

Finally, we give our main result on minimum level spacing. Let {L^a}a=i,..., 2 |A| denote 
the eigenvalues of H. 


Theorem 1.3. (Theorem 3.9). There exists a constant b > 0 such that for'y sufficiently 
small, 

p( mm\Ea — E/sl < S) < |A| exp (— 6 | logyl^^'^l log5|^/^) , (1.7) 

\ (y.^13 / 

for any rectangle |A| and any 0 < h < 

As this is a result about microscopic level statistics, we are interested in the behavior 


as 5 —)■ 0. This bound tends to zero more slowly than (1.5), but still faster than any 
power of 1/1 log(5|. 


1.4 Fundamental Lemma 


We state here a lemma on Schur complements that will be used often in our analysis. 


Lemma 1.4. (Fundamental Lemma) Let K be a {p + q) 


block form, K = 

Assume that ||(il — E)~^ 
with respect to \: 


A B 
C D 


, with A a p X p matrix, D 
< e-\\\B\\ < 7 , lie'll < 7. 


Ex = A- B{D - X)-^C. 

Let e and ■j/e be small, and |A — i?| < e/2. Then 


X {p + q) symmetric matrix in 
a q X q matrix, and C = B^. 
Define the Schur complement 

( 1 . 8 ) 


1. If ip is an eigenvector for E\ with eigenvalue \, then — \)~^C(p) is an 

eigenvector for K with eigenvalue X, and all eigenvectors of K with eigenvalue X 
are of this form. 

2. The spectrum of K in [E — e/2, E + e/2] is in close agreement with that of Ee 
in the following sense. If Xi < X 2 <...< Xm are the eigenvalues of K in 
[E — e/2,E -|- e/2], then there are corresponding eigenvalues Xi < X 2 < ... < Xm 
of Ee, and |Ai - A/ < 2(7/e)^|Ai - E\. 


Note that the degree of closeness between the two sets of eigenvalues improves the 
closer one gets to X = E. 

Our goal, then, is to iterate the fundamental lemma, using appropriately chosen 
subspaces, so that the spectral window width e tends to zero, with 7 /e tending to zero 
as well. Eventually, the width will be so small that at most one eigenvalue will be 
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inside. At that point, the decay of the eigenfunction will become manifest, as will the 
separation of the eigenvalue from the rest of the spectrum. 

The operator —{D — X)~^C plays the role of expanding a solution of a low dimen¬ 
sional space to a higher dimensional space. Such factors will be composed with each 
other to obtain the proper eigenfunction on the full Hilbert space; hence we call them 
eigenfunction-generating kernels. 

Proof. The hrst statement is a standard fact about the Schur complement. If ip is 
an eigenvector for F\ with eigenvalue A, then 

/ ^ \^f{A-X)p-B{D-X)-^Cp\ 

^ \-{D-X)-^CpJ \ Cp - {D - X){D - X)-^Cp J 

= ( ) = 0. (1.9) 

Conversely, if {K — X){p, (p) = 0, then Cp + {D — X)p = 0, which has unique solution 
p = —{D — X)~^Cp. Substituting in, we obtain {Fx — X)p = 0. Note that this implies 
agreement between the multiplicities of the A-eigenspaces for H and Fx. 

For the second statement, write 

Fe-Fx = -B{D - E)-^C + B{D - X)-^C = B{D - E)-\X - E){D - Xy^C. (1.10) 

Noting that |A — i?| < e/2 and dist(spec iA, i?) > e, we have that ||(iA — A)“^|| < 2/e, 
and thus ^ 

\\Fe-Fx\\ < \\X-E\. (1.11) 

By Weyl’s inequality, the eigenvalues of Ex and Fe differ by no more than 2{'y/e)‘^\Xi—E\ 
when shifting from Ex to Fe = Ex + {Fe — Ex). □ 

2 Fixed Energy Procedure 

2.1 First Step 

Let us investigate the spectrum of El in the vicinity of E. The energy E can be any 
real number (although since H is bounded, there will be no spectrum near iA if |iA| is 
large). For the time being, E will be treated as a hxed parameter. We emphasize that 
we are deriving probabilistic bounds on resonances by way of the geometric/algebraic 
behavior of the eigenvalues with respect to the random parameters. Let us observe this 
in the simplest possible case. 

It will be helpful to introduce some terminology. The rectangle A is a collection of 
positions in labeling the coordinates in the hrst factor in "H = (8)C^. Coordinates 

for the full PL will be called sites] thus there are |A| positions and 2|A| sites. Let A* 
denote the set of sites. More generally, let X* denote the set of sites at positions 
X E X C A. It can be visualized as two copies of A in two layers: two sites at each 
position. 
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The eigenvalues of are ^.t^i where 


U = ^ul + 1 . ( 2 . 1 ) 

We say that x is ^-resonant with E if either eigenvalue is within e of E. Sites associated 
with non-resonant positions are also said to be non-resonant. The eigenvalues are split 
by at least 1, hence for e: small only one can be close to E. A direct calculation shows 
that 


P{U e [\E\ - e, \E\ + e]) < P{t^ G [1,1 + 2e]) = 2Ve + < 3 ^, for e < 1, (2.2) 

assuming a uniform distribution on [-1,1] for (The maximum is achieved when 
|£'| = 1 -|- e.) Thus Si/i is a bound for the probability that x is an e-resonant position 
with respect to E. We take e = 7 "^ with 0 = |. 

Let denote the set of e-resonant positions x G A. It can be broken into compo¬ 
nents, where connectedness is based on nearest-neighbor connections. Let us estimate 
Pi]l, the probability that x and y are in the same component. For positions x and y 
to be connected, there must be a self-avoiding walk from x to y, consisting of nearest- 


neighbor steps between resonant sites. By (2.2), the probability that all the sites in a 
given walk are resonant is bounded by {3y/e)^, where m is the number of positions in 
the walk, including endpoints. If we add a combinatoric factor 2d per step, we may 
convert the sum over walks with m positions into a supremumj^ A combinatoric factor 
2™ allows us to £x m. We obtain 


PW < (i 2 dV^)- < {12d^y^-y\+\ (2.3) 


as m — 1 must be at least \x — y\, the minimum number of steps to get from x to y. 

Block form of the Hamiltonian. Recall that is the set of resonant positions. 
Putting = A \ R^^\ we have associated sets of sites R''^^*, These index sets 

determine the block form of the Hamiltonian: 


H 


71 ( 1 ) ^( 1 ) \ 

( 7 ( 1 ) /)(!) j ’ 


(2.4) 


with representing the restriction of H to R^^'^*, and representing the restriction 
to (or, to be precise, these are the restrictions of H to the subspaces corresponding 

to those index sets). Note that contain off-diagonals only between sites at the 

same position or at adjacent positions; positions which belong to different components 
in i?d) clo not interact. 

We may write 

L)P) = IRW+ (2.5) 

with block diagonal (with blocks at each position x) and with equal to 
the projection of V to i?d)c pet Vxy'^ denote the block matrix element of between 

1 Combinatoric factors are a convenient bookkeeping device when estimating sums. If ^ 

for some positive constants Cq,, then < sup^, |Xq,|cq. We call Cq a combinatoric factor. 
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positions x and y. For x,y E it is equal to 'jJxyh, see (1.2). Thus it is nonzero 


only for adjacent positions x and y in in which case it equals 7 / 2 . 

Convention. We use ordinary absolute value signs | ■ | to denote the norm of a 2 x 2 
matrix such as the block matrix element Vxl'^. This is to distinguish it from the usual 
notation || ■ || for matrices indexed by A* or some subset of A*. Using this notation, we 
have that |14y^| = 7 for x^y adjacent positions in 


Random walk expansion. We would like to reduce (2.4) by Schur complementation: 


ird) = 7i(i) _ 


( 2 . 6 ) 


Let us develop a random walk expansion for — A)“^. 

By construction, dist(spec U) > e. Let us assume that |A — U| < e/2, so that 
||(pp(i) _ A)“^|| < 2/e. The off-diagonal matrix has norm no greater than 2^7 
in dimension d (it is a symmetric matrix with row sum equal to 2 ^ 7 ). Therefore, for 
s = ■y't = sufficiently small, — A)“^ exists and can be expanded in a Neumann 
series: 


(D(^) - A)-^ = (lU^^) - A)-^ + (lU^^) - A)-V(^)(1U(^) - A)-^ + .... (2.7) 


This leads to a random walk expansion with steps determined by block matrix elements 
of 


m—1 


= E n w 


Xj^l * 


( 2 . 8 ) 


gi'.x^y i=l 


i=i 


Here gi = {x = X\,X 2 , ■ ■ ■ ,Xm = y} is a random walk with each Xi, Xj+i adjacent 
positions; return visits are allowed. For such a walk we write its length as \gi\ = 
m. Note that each factor on the right-hand side of (2.8) is a 2 x 2 block matrix 
element. As each eigenvalue of — A is at least e/2 in magnit ude , we have that 
|[(pp(i) _ X)~^]x^x^\ < ‘2/e for every Xi E gi. Thus a term in the sum (2.8) with l^fil = m 
is bounded by 7 ™'(e/ 2 )“^™+^b We may assign combinatoric factors 2™ to fix the length 
of the walk and 2d for each sum over Xj] altogether the combinatoric factor is 2"^{2d)"^~^. 
Thus we obtain a bound 


[(D<‘l - A)-‘] 


xy 


< 


4 


\ e / 


(2.9) 


as m — 1 must be at least \x — y\, the number of steps needed to walk from x to y. Here, 
and in what follow s, Cd deno tes a constant that depends only on the dimension d. 

Observe from (2.7) and (2.8) that the operator (D^) _ yj-i connects all sites in 
7 ?b)c* Therefore, we see that resonant sites ^re conn ected through the chain of 
operators in the second term of (2.6). Note that and 

contain off-diagonal interactions from V. Thus the walks in (2.8) are extended with 
two additional steps VJJL 


rjXj_f.i- 


As in the arguments for (2.9), we obtain 




xy 


,(£n) 


\x-y\y2 


■ e. 


( 2 . 10 ) 















In a similar manner, we can estimate the eigenfunction-generating kernel for the hrst 
step: 




xy 


< 


(f) 


h-?;|vi 


( 2 . 11 ) 


The graphical expansion (|2.8[) indnces similar expansions for these matrix elements. 


2.2 Resonant Blocks 

In the step we use a length scale Lk = to determine connectivity. That is, 
in the hrst step, Li = 1 and two positions are in the same component if they are 
nearest neighbors. In steps /c > 2, we inherit a set of resonant positions and two 

positions are in the same component if they are within a distance Lk- Thus, components 
at scale Lk are subsets of that can be spanned by connections of length up to 

Lk- A component of cardinality n has 2n sites; each of these is effectively connected 
to all the others of the component. Distance will always be measured in terms of the 
original metric |a:| = 1^*1 

Definition 2.1. Let B be a component of on scale k with k > 2. We say that 

B is isolated on scale k if 

1. its volume (= n{B) = number of positions) is no greater than , 

2. its distance from every other component on scale k is at least ALk, 

3. its diameter is < Lk- 

Here ^ is an exponent governing fractional exponential decay of graphical bounds. 
From the previous step, we have the Schur complement matrix 

( 2 . 12 ) 


with bounds similar to (2.10) that are proven below in Theorem 2.7 (and that can be 
assumed by induction). It operates on the space indexed by elements of Each 

of the matrices in (2.12) is given by a random walk expansion as in (2.8) in terms of 
the corresponding matrices from step k — 2. See Section 2.4 for details on the random 
walk expansion; for now, we need only some of its general properties. 

In each step we take the Schur complement of the current renormalized Hamiltonian. 
In the hrst step, is a renormalization of the original Hamiltonian H] the interactions 

of the hrst step generate walks, or graphs constructed of connections in the original 
Hamiltonian. Further renormalizations ... involve walks of the preceding 

Hamiltonian’s walks; we will refer to these nested walk structures as multigraphs. Thus 
the Schur complement F^^ is dehned inductively in terms of (/c —l)-level multigraphs, 
each being a walk with steps given by {k — 2)-level multigraphs, and so on down to the 
hrst level. By construction, these graphs sample the potential only for positions 
X that the hrst-level walk passes through; thus there is dependence on only for 
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X ^ Consequently, the dependence on u^, x G R{k-l) 

is only through the 

leading term in specihcally the matrix projected to (Recall 

that y 4 d) is the original Hamiltonian H projected to In particular , th e values 

Ux G do not appear in any of the block inverse operators used - see ( 2 . 8 ) and its 


step version (2.55) 


One may visualize the last term in ( 2 . 12 ) as a set of graphs that exit R^^ (the 
matrix and then wander in ]^{k "i) ^ R{k 1 ) matrix b _ yj ij before 

returning to (the matrix Thus we see that agrees with H on 

R{k-1)* 

up to small terms from the random walk expansions, all independent oi Ux, x ^ 
Rik-i)_ Hence is equal to a direct sum of H projected to the R^^ components 

plus small 0{'y^/e) corrections (independent of u^, x G That is, we may write 


— (Bshs + Sfc-i, 


(2.13) 


where hs is the projection of H to the com ponen t B in a nd Sfc_i = /e) 

is independent of u^i x G For k = 2 (2.13) holds by (2.10) and indeed we will 

generalize it for fc > 2 in the induction procedure - see Corollary |2.8 Furthermore, 
there is no off-diagonal dependence on Ux, x G 

We wish to show that isolated components are unlikely to contain spectrum near E. 
We could work with projected to an isolated component B of but in order 

to preserve independence we truncate the random walk expansions in including 

only multigraphs that remain within a distance < Lk -2 of B. Multigraphs neglected 
therefore travel a distance > 2Lk-2 = L^-i after returning to B or reaching another 
block. Let us designate such multigraphs nonlocal at scale Lk-i, they visit points at a 
distance > Lk _2 from B. Let us write F^~^^ for the renormalized Hamiltonian with 
all nonlocal multigraphs dropped. The truncated Hamiltonian is block diagonal due to 
the separation condition: F^ 


(fc-i) _ 


= ®bF^ ^\b)- As above, we may write f]^ ^\b) = 


liB -l- Sfc_i(F) for the truncated f|*^ projected to block F, and Ek_i{B) = 0{'y‘^/e). 


We prove decay at rate for graphs of length L, with 'tp = \ - see Theorem 


2.7 


below. We prove below a bound on the probability of spectrum within £ of E. As 
this window is much larger than the size of terms neglected (recall that £ = 7 '^ with 

1 r 

0 = ^), movement of the spectrum due to terms of order 7 is insignihcant. 

Definition 2.2. Let B be a component of on scale k with k > 2. We say that 

B is resonant on scale k if it is isolated on scale k and z/dist(spec F^^~^^(F), F) < 
Ek = E k, Here we refer to the truncation Ff, —>■ ©bF^ [B), which results from 
the deletion of all multigraphs that extend a distance Lk -2 or farther from any block B. 


Following a similar plan as in the hrst step, we dehne the new resonant set R^^"^ by 
removing from all the non-resonant isolated components. Thus 

R{k) ^ ^(fc-1) y y ^2.u) 

a: Bq, is isolated but not resonant on scale k 
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Proposition 2.3. Let B be an isolated component of on scale k with k >2. The 

.3^/5 yo 

probability that B is resonant on scale k is less than . 

Proof Let us enumerate the positions in B as {xi,..., x„}. The block B is connected 
on scale Lk-i (he. using connections between Xj and Xj of length up to Lk-i. Let us 
consider the 2n x 2n matrix F^~^\b) as a function of u = {ux ^-,..., with all other 
m’s held hxed, which means that Sfc_i(i?) is hxed as well. Thus, as explained above, 
F^~^\b) agrees with F[ up to u-independent terms of size 0 ( 7 ^/e). We consider here 
X = E, since by the fundamental lemma, control of the spectrum of F^ ^\b) allows 
us to control F^^ ^\b) for A near E. Recall that H has Ui or —Ui on the diagonal, 
and it has off-diagonal matrix elements connecting adjacent positions and connecting 
the two sites at each position. The off-diagonal matrix entries are either 7 , 1, or 0. It 
should be clear, then, that A(u) = det(T’g^~^^(i?) — is a polynomial in u of degree 
2n. We may therefore apply standard results bounding the measure of sets where the 
polynomial is small. 

Noting that each Ux G [—1,1], and since tx = -|- 1, the spectrum of H is clearly 

limited to the range [—\/2, a/2] - 1 - 0 ( 7 ). Therefore, we can assume G [—|, |]; this will 
cover all of the spectrum of H. 

With E G [— 1,|], we need to establish in a quantitative sense that A(u) = 

det(T’g^~^^(R) — E) is not the zero polynomial. Thus we look for a value of u for 

which A(u) is bounded away from 0. Take Uq = (2,..., 2). Then each tx^ = a/5, and so 

{tx - E){-tx - E) = - 5 < I - 5 = - A. (2.15) 


We must allow for 0 ( 7 ) movement of the eigenvalues due to the terms other than 
®ihx^ in F^~^\b) (here we use also the decay away from the diagonal - see Theorem 


2.7 below - and the fact that the norm of the perturbing matrix is bounded by the 
maximum absolute row sum). Still, we maintain a lower bound |A(uo)| > 2 ”, after 
taking the product over the n 2 x 2 blocks and allowing for the small changes in the 
eigenvalues. 

We use the Brudnyi-Ganzburg inequality, in the form stated in [B]: 


sup IpI < 
u 



K 

sup \p\. 

U) 


( 2 . 16 ) 


Here f/ is a bounded convex subset of M”, p is a polynomial of degree at most k and u 
is a measurable subset of U. In our case, we let U = [—2,2]"', k = 2n, and 

(u = {u : |A(u)| < with Sk = ■ (2-17) 

We have shown that sup^/ |A(u)| > 2"' and so 

^ ■ An\U\ < ■ nA^. (2.18) 
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Since is distributed uniformly on [-1,1], this leads to bounds 
F(|A(u)| < . 9 n 2 ” 


< e 


^3V>/5 


/2 . 


2^/5 lo 


(2.19) 


where in the second line we have used the bound n < from Definition 

well as the smallness of e = Note that both E and the eigenvalues of 


2.1 

(n) 


''k 

are m l —|, |], and therefore the eigenvalues of — E are no greater than 3 in 


(as 

{B) 


magnitude. Hence (2.19) implies the statement of the proposition. 


□ 


2.3 Probability Estimates 


We need to work on a multiscale version of the percolation estimate (2.3). As we 


investigate the percolation of resonant blocks, we face a situation in which blocks of 
flW ] 2 iay inherit probability bounds from step fc — 1 or they may receive a probability 


bound from Proposition 2A Thus we make an inductive dehnition that keeps track of 
the probability bounds available for blocks at any scale: 


P^^\B) e 


at/'/s . 


s^k -2 jf jg isolated on scale k — 2] 

nm'ip(A:-i)^^.)^ Otherwise. 


( 2 . 20 ) 


We emphasize that the P are not themselves probabilities but stand in for upper bounds 
on the probabilities of resonant clusters. Here n{B) is the number of positions in B, 
and Bi,... ,Bm' are the subcomponents of B on scale k — 1. Note that each factor 

. 3 - 0/5 

3y£ or £ fe-2' (corresponding to a bound on the probability that a site or block is 
resonant) carries forward to the next scale, until a scale k is reached with B isolated 
on scale k — 2. In that case, all the factors from subblocks of B from previous scales 


are replaced with the bound from Proposition 2.3 This corresponds to replacement of 
bounds by a better one from Proposition 2.3, when it is available (he. if a component 
is isolated on scale k — 2). 

Next let us define a weighted sum of the probability bounds P^’^pB): 




B containing x and y 


P^^\B)e 


-qkn{B) 


( 2 . 21 ) 


Here {qk} is a decreasing sequence of numbers with a strictly positive floor q > 0; 
it will be specihed in the next theorem. We prove the following “energy-entropy” 
bound, similar in spirit to ng. It is a purely geometric-combinatoric result about the 
components B that are generated by our dehnitions. 


Theorem 2.4. Given '0=3, let x = V'/S- Then for s sufficiently small, and qi = 


92 = |. 93 = h, and qt = 9i_i(l - L^*!p) for k>i, we have 


Q'tl< 


9fe(h-y|vLfe_2)^ 


( 2 . 22 ) 
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Proof. To treat the case /c = 1, we need to extend the proof of (2.3) to incorporate 


decay with n = n{B), the volume of the component containing x, y. The walk from x to 
y can be extended to a branching random walk (he. a spanni ng tr ee). The number of 
such trees is bounded by cJJ. The weighting in (2.21) is counterbalanced 

by our bound on the probability that n positions are resonant, {3y/e)"‘. Noting that 
n > |x — |/| + 1, we may extract the decay in (2.22), with a decrease of an additional 


1/5 in the power of e. This leaves a bound of . The sum on n is less than 1, 


and hence (2.22) holds for k = 1. 


Let us now assume (2.22) for A; — 1. Let us write 


f){k) ^ , f)ik),2 


'x,y 




(2.23) 


with the hrst term giving the sum over blocks B that are not isolated on scale k — 2, 
and the second term giving the sum over the ones that are isolated. We consider 
hrst and note that P^^\B) is given by the second alternative in (2.20). Recall 


that R is a component of and that it is formed by joining together components 

of with links of length < L^. Hence there must be a sequence of positions 

X = xi,yi,X2,y2, ■ ■ ■ ,Xm,ym = 2/ such that each Xi^yt lie in the same component Rj 
of and | Xi — yi-i\ < Lfc. The remaining subcomponents of R may be denoted 

Rm+i, • • •, Rm'- Connectivity on scale k implies that there must be a tree graph T 
connecting all the components of R. Each link of the tree graph is a pair (xi,|/j_i) as 
above with 2 < i < m or a pair {xj,yj) with m + 1 < j < m', Xj G Rj, ijj G B^-^j), 
and r(j) < j. Thus we choose a tree graph that extends the unbranched connection 
from Ri to Bm- (See Figure [^) Again, because of connectivity, we can assume that 

\xj yj\ ^ Lfc. 

We may write 


m 

Q'P' <T.I1 E 11 


-qkn{Bi) 


(2.24) 


m,m' T i=l 

consistent with T 


where we have written n = of volumes of the Bi. With 

some overcounting, we may proceed through the tree graph, summing over vertices 
Xi,yi,Xj,yj and components, requiring only that each Rj contain a specihed set of 1 
or 2 vertices. In this way, we will be able to apply the inductive hypothesis from step 
k-1. 

The tree graph T may be broken up into the chain connecting Bi,..., Bm plus 
individual tree graphs Ti,..., Tm with roots at the corresponding blocks. We now focus 
on a single Bi (which we take as hxed) and sum over R and the associated blocks (other 
than Bi) that are linked by Tj. Let us dehne for an arbitrary block R on scale k — 1\ 




n 


T: depth(T)<p 


Bi 


i=l 


(2.25) 


consistent with T 
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B\ B2 -B3 



Bi 


Figure 1: A tree graph of blocks connecting x and y. Blocks form the spine 

and Bi is a leaf. 


Here depth(T) is the largest number of links in T that are needed to reach a vertex, 
starting at the root. Then we can write 


m=l yi,X2,y2,---,ym-l,Xm i=l LSi containing 


(2.26) 


Here the inequality results from the overcounting associated with treating each of the 
Ti as independent of the others. Let us put 

£ = k>2. (2.27) 

Lemma 2.5. With X = § o,nd e sufficiently small, 

K^p\B) < (2.28) 

Proof. We use induction on p. When p = 0, we have the tree with no links, so ( |2.28[ ) 
amounts to the statement 1 < For larger p, we may use the recursion 


00 r 

< En 


r=0 i=0 


- Bi linked to B 


(2.29) 


Again, this is an upper bound, as we are treating each of the r branches of T as 
independent of the others. There is also overcounting due to permutations of the 
branches, and this has been compensated by the factor 1/r!. 
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The sum over Bi linked to B may be bounded as follows. By connectivity, there 
must be a pair of points x, y with y E B and x E Bi and \x — y\ < L^. The sum over 
y leads to a factor n{B). There are no more than {2Lk + choices for x, given y. 
Furthermore, we may bound 

(' 2 . 30 ) 


using (2.28), since 


-{qk - qk-i)\loge\ + f > 

> 0 


0 , 


loge\qk-iLi^^l^\ k>A 
k = 2,3 


(2.31) 


for e small. (This follows from the expressions for qk in Theorem 2.4, in particular 
the relation qk-i — qk = qk-iLk-^^ ^ > 4; the increase from -0/20 to -0/15 is to 
accommodate the e term on the left-hand side.) Inserting (2.30) into (2.29), we obtain 

K^P\B) < exp {{2Lk + l)"g^,"'^n(5)^ 

< exp [{2Lk + l)‘^£'''=-i^^-3n(.B)) < 


(2.32) 


after applying the induction hypothesis (2.22) to ■ This closes the induction 

and completes the proof of the lemma. This “polymer expansion” construction is a 
standard way to treat tree graph sums when the vertices have spatial extent, see for 
example [18] . □ 

After using the lemma to sum over all the branches of the tree, we are left with the 
“spine,” that is, the blocks i?i,... ,i?m that lead from x to y. From (2.26), we obtain 
the following bounds: 


Qih" < E E n 


m=l yi,X2,y2,---,ym-l,x:m *=1 LSi containing 


p{k-l) 


(2.33) 


< 


Bl containing x and j; m >2 yi,X2,y2,---,ym-l,^m i=l 


The effect of the tree graphs has been subsumed into factors by Lemma 2.5 


The hrst term of the hnal expression bounds the case m = 1. For the terms m > 2, we 
have used the fact that 


Bi containing Xi^yi 


(2.34) 


which follows from (2.21),(2.31). 
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Consider the term m = 2, and let e = £'^'=- 1 . We need to be carefnl with the snms 
on ?/i,X 2 in order to preserve the form of the decay estimate (2.22). With x,y hxed, 
note that 

min{|x -yi\,\x 2 - < CdL'^. (2.35) 

y\,^2 

This can be seen by considering two cases, depending on whether the minimnm is \x—yi\ 
or \x 2 — y\-i and noting that \yi — X 2 \ < L^. As a conseqnence, a combinatoric factor 


CrfL^min{|a: -yil\x2- 


(2.36) 


snffices to control the snms. Here Lmin = min{La,Lfe}, Lmax = niax{La,Lfe}, with 
La = \x — yi\\J Lk- 3 , Lb = \x 2 — y\y Lk- 3 . Hence we may bonnd the m = 2 term by 

E < ^up £«£«, ( 2 , 37 ) 


271 ,2:2 


yi,^2 


yi,x2 

Note that 

1^ 2/| N L]^_2 E La T T/j T Lb T L^: T L^^yii (2.38) 

since as indicated above, \x 2 — yi\ < by connectivity of B. Let us write 

CdL^^+^e^^e^i = ^2.39) 

which holds for small e. We claim that 

t-^min + -^max > (-^min + -hfc + Lmax)^- (2.40) 

The worst case for this inequality is when La, Lb are at their minimum possible values, 
Lk -3 (as one can check by differentiating). Then it reduces to the inequality 1.8 > 10^, 

or y < log^g 1.8 ~ .255. Recalling that '^ = | and y = 'ip/3, we can conhrm this 

inequality, and then (2.38),(2.40) imply a bound 


1 ~(|a;-y|VLfe_2)^ _ 1 (jfe_i(|a;-y|VLfe_2)^ ^ 1 (jfe(|a;-y|VLfe_2)^ 

4^ 4^ ^ 4^ 


(2.41) 


on the term m = 2 . 

The terms m > 2 may be handled inductively, since the bound just proven is suf- 

hcient to reproduce the argument when summing each pair yi,Xi+i in turn, and we 

obtain a bound on the term. 

Now consider the term m = 1, i.e. the hrst term in (2.33). This corresponds to a 

situation where x,y are already contained in the same block on scale k — 1 . As with 

ik—V\ 

the m >2 terms, we have a bound given by Qx,y on that scale: 




(fc-i) < 

x,y — ^ 


qk-l{\x-y\y Lk-zY 


(2.42) 


We need to improve this by replacing Lk -3 with Lfc- 2 - For this, we can assume that 

\x-y\< Lk- 2 - 
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First, suppose the scale k — 1 block containing x,y is composed of 2 or more scale 
k — 2 blocks. We obtain a bound which is more than sufficient, because we receive two 

(Here we 
will allow 


or more factors of from (2.22), which leads to 


use the fact that X = S) so 2,L\_. = 2^ and the excess over Lj 


us to control the sum over relative positions of the two blocks as in the arguments for 
(2.35)-(2.41) above.) We should also consider the case k = 2, for which we hav e at least 


a factor for a single block on scale 1, as explained in the proof of (2.22) for 

k = 1. 

On the other hand, suppose the {k — l)-block containing x, y consists of a single 
{k — 2)-block. If A; = 3, then A: — 2 = 1 and , so the desired bound 

holds. If not, then at least one of two cases must hold: n > i diam(i?) > Lk_ 2 . 
(Since the block does not join up with other blocks after two jumps in scale Lk _2 —)■ L^, 
its distance from other blocks must be at least 4Lfc_2. As we are working on the case 
where B i s no t isolated on scale k — 2, we must violate one of the othe r con ditions of 
Dehnition 


2.1 


provided k — 2>2.) For the case n > have by (2.31) a bound 


p(fc-i)(5^)^-9fe-ir^(Si)exp (^-|fog£|gfc_iL":^’/^®n(Hi)) (2.43) 

Bi containing x and y 


on the hrst term of (2.33). The sum over Bi is bounded by (2.42), and we pick up an 
extra factor 

Y — i/^/lS 

fr1k-l^k_2 ^ 


’li < g1k-i^k-2 = (2.44) 

which includes the desired improvement Lk _2 —)■ Lk_ 2 . Finally, consider the case of a 
block on scale k — 2 with diameter > Lk- 2 - If this is the case, we must be able to hnd 


x',y' in the block with \x' — y'\ > Lk- 2 - Then we can use (2.22), which becomes 

/p(fc-2) ^ < f<lk-2Ll_2 

^ x',y' — ’ 


(2.45) 


and after summing over < { 2 Lk -2 + 1)“^ possibilities for x' and for y', given that x, y are 
in the block, we obtain a bound as in the other cases. (The drop qk -2 — >■ qu-i 


allows us to use the additional factor ks k -2 fo control the sums - see (2.46) 


below.) Altogether, we have considered up to three cases for any particular k. Let us 
allow \x — y\ > Lk -2 again. We obtain the following bound on the m = 1 term: 


3^9*;-l(h-y|VLfc_2)^ < ^^qki\x-y\yLk- 2 )^^qk-iL^.tl!(^°LI _2 ^ l^qk{\x-y\\/Lk- 2 )^ 


(2.46) 


(k') 2 

We return to a consideration of Qx,y , which consists of terms from (2.21) that are 


isolated on scale A: — 2, which means that A: > 4. In this case we have from (2.20) that 


P^^\B) = We need to sum over all possible blocks of diameter < Lk -2 and 

weight by to obtain Q^x}y‘^. We obtain a bound 


£^^"-2 p{2Lk-i + l)'^^e:“9fc»^ ^^7^2 B < \£^kLl_2 = i^qk{\x-y\'-JLk-2Y^ 


rSV'/S 


(2.47) 


using dnlog( 2 Lfc _2 + 1) < CduL'^^.^ < CdL'^j^J^ and qkU < qkL^^-^ < ^^-2 
qk < ^ for A: > 4, which can readily be checked from the recursion). Note that we are 


3y/5 


-2ip/5 


-sy/s 
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using the Proposition |2.3| bound to control the sum over all substructure of B. This is 
because Proposition |2.3| cannot be used simultaneously for a block and for its subblocks, 
due to lack of independence. However, the bound for the whole block is much better 
than what had been obtained earlier for subblocks, and (2.47) demonstrates that it is 
adequate for our purposes here. 

To complete the proof of the theorem, we gather the bounds (2.46), (2.47) and the 
bound proven above on the term, m > 2. The desired bound 

(2.22) follows, since | + | + Y.m >2 2”™ = 1- D 

We may now estimate the probability that x and y belong to the same resonant 
block on scale k. 


Corollary 2.6. Given 'tp 
g > 0 such that 


X = ip/3, let e he sufficiently small. Then there is a 


ffik) ^ pqi\x-y\VL^,_2)>^ 

’ x,y — ^ 


(2.48) 


Proof. By construction, the probabilities and in (2.20) correspond 

to independent events, because whenever a block is resonant on a scale fc, its probability 
bonnd from Proposition 2^ replaces all previous probab ility factors for its subblocks. 


j Sip jh I I 

So any two snbblocks with a factor eG m appearing in (2.20) are separated, ensnring 
independence. (Independence is a result of the truncation of the random walk expansion 
to form including only graphs that remain within a distance < Ta :-2 of B, 

so that there is no overlap between sets of variables u on which each resonance event 
depends.) At each length scale, isolated blocks must be resonant (otherwise they wonld 


/3 Proposition 


be removed before forming blocks at the next sea' 
a probability factor e 


2.3 


) and hence can be associated with 
Thus all the probability factors in 
(2.20) correspond to independent events that must hold if H is a component of 
Hence P^^\B) is a bonnd on the probability that H is a component of 

Since the weights in (2.21) are bonnded below by 1, we see that (2.22) becomes 
a bound on the snm of PP\B) over B containing x,y. The corollary follows, since 

gfc > g > 0. □ 

The underlying mechanisms at work in the proof of Theorem 2.4 are as follows. 
First, since we are only seeking fractional exponential decay, the procedure can tolerate 
giving up some fraction of the decay distance with each new length scale (which happens 
because of the gaps that were required between resonant components). Smallness of the 
probability of a component comes in one of two ways. If a component is isolated, with 
small volume and diameter, the basic proposition on resonance probability can be used 
to bring in new smallness to control its greater positional freedom - it gives an estimate 
that decays as a fractional exponential of the length scale. On the other hand, if the 
component’s volnme or diameter is large, or it is not too far from other components, 
the decay can be obtained from previous scales. We need to be careful not to give np 
too mnch volnme decay for this purpose in each step, so that we can obtain bounds 
that are nniform in the step index k. 
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2.4 Random Walk Expansion 

The random walk expansion is based on the Schnr complement at level k. We pnt 


^ ( /K'') BK') \ 

I (j{k) jj{k) \ 1 


(2.49) 


where the blocks are determined by the decomposition of into (npper-left 

block) and \ R^^'^ (lower-right block). We may now dehne the k}^ Schnr comple¬ 

ment for |A — i?| < ek/‘2'. 


p{k) ^ ^{k) _ ^(fc)p(fc) _ ( 2 . 50 ) 


Let ns write 

^ pp(fc) (2.51) 

where is block diagonal, each block being ^^(-B) for some B. This means that 
consists of the long graphs not inclnded in F^~^\b)^ prodncing matrix elements 
both within blocks and between blocks. 

We show below in Theorem 12.91 that 

||Ff-b( 5 ) _ F^^-^\b)\\ < Cd-^\\- El (2.52) 

which is less than £fc/ 6 , becanse |A — < 6^/2. Since all the blocks of \ R^^'> are 

non-resonant, 

dist (^spec F^~^\b), >ek = , (2.53) 

and so 

||(Vp{fc) _ A)-i|| < (2.54) 

Hence, as in the hrst step, we may expand — A)“^ in a Nenmann series: 

m m—1 

- A)-‘w = ^ (2-55) 

Qk-x^y *=i i=i 

Here gk = {x = xi, aii, X 2 , X 2 , ■ ■ ■, Xm, Xm = y}, with each a;*, Xi in the same block. Note 
that Vx^y is given by a snm of graphs contribnting to F^ — ®bF^ where 

(2.56) 


Thns, when it is nsefnl, we may nnwrap any matrix element into a snm of mnltigraphs, 
since each step at any level consists of a s um of graphs on the previous level. Thus, the 
matrices in (2.56) are themselves blocks of F^~‘^\ which 

is given by We may continue down to graphs at 

level 1. 
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We will also prove bounds on the matrices that generate the eigenfunctions. Recall 
from the fundamental lemma on the Schur complement that if is an eigenvector of 


with eigenvalue A, then 




(2.57) 


ik—V) 

is an eigenvector of F) with the same eigenvalue. This process may be repeated to 
extend the eigenvector all the way out to the original lattice A*, that is, to produce 
an eigenvector of H. Let us write 


^(»l = G'‘V 


(k) 


and then we may give a multigraph expansion for G 


(k) 


(2.58) 

in the same manner as was 


just c 


escribed for Indeed, the same operators appear when unwrapping 


(2.57). The matrix can be thought of as half of F'^’ since it descends to the original 
lattice, but does not climb back up to 

One may visualize the multigraph expansion on the original lattice as an ordinary 
random walk of nearest-neighbor steps, except that when it enters a block on level A;, 
there is a matrix element of — A)“^ that produces an intra-block jump. This 

structure is present already at the hrst level, see ( 2 . 8 ); but there we had only single¬ 
position blocks. The larger blocks that appear in later steps create gaps in decay, since 
there is no decay from Xj to x,- in the following estimate for a walk passing through a 
block B of 


(fc) 






(2.59) 


Xi 


Here we use (2.54), and note that there are just n{B) choices for x,, given Xj. As 
the block is isolated on scale k, its volume is bounded by ■ The main issue with 
controlling the multigraph expansion is making sure that the gaps in decay and the 
large factors from the bound (2.59) do not spoil too much the exponential decay proven 
in step 1. Our constructions ensure that there are gaps of size 4Lfc between blocks at 
scale k. Hence, whatever rate of decay is proven at scale k — 1, some small fraction will 
be lost at scale k. This is why we are led to fractional exponential decay estimates. 
Keep in mind that as far as the step k random walk expansions are concerned, the 
conhguration of blocks up to that level R^‘^\ ..., R^^'>) is hxed. 

We now state our main theorem on graphical bounds. Let Sx,z,y denote the sum of 
the absolute values of all multigraphs for that go from x to y and 

that contain z. Here, x^y are positions in R^^\ and z is a position in R^^'^ \ R^^ We 
say that a multigraph contains z if any of the blocks that it passes through contain z. 

Theorem 2.7. With -0 = |, 0 = and e = sufficiently small, assume that |A —F| < 
Ek/^. Then 

< r^[i\x-z\+\z-y\)y lW . 2“^ (2.60) 

< ry[(h-^l+h-y|)vi]’^_ (2.61) 

j<k 


E 


^(fc) 

x,z,y 

SU) 

^x^z,y 
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Proof. For the case k = 1 we obtain 


sp) ^/^yi^-"i+i"-yi)v2 _ 


x,z,y — \ g J 




(2.62) 


as in th e pro of of (2.10). Let L = {\x — z\ + \z — y\)\J 1, then one has a power |(L V2) + ^ 
of 7 in (2.62). This is greater than the required power for (2.60), as one can readily 


verify, checking the hrst few values of L by hand. Hence for 7 sufficiently small, we 


obtain (2.60) 


Let us now assume (2.60) for k — 1. Using (2.54), one may expand _ 


y)-i(7(fc) to obtain a graph connecting x to y traversing m inte rmed iate blocks. In the 
hgure, each block has a matrix element of — A)“^ as in (2.59). The extra point 

z must lie in one of the sections of the graph (as in Figure]^ or in one of the blocks. 
Because the blocks are isolated, and x, y are in each section of the graph traverses 





Figure 2: A graph from x to y in the step k random walk expansion 


a distance > and the blocks themselves have diameter < L^. We sum over the 
pairs (xj, Xi) in turn in a manner similar to what was done for the probability estimate. 
In Figure]^ for ex ample , the sum over x, may be controlled with a combinatoric factor 
Cd\x — Xi\^^^ and (2.59) controls the sum over xi and results in a factor . 

Applying the inductive hypothesis to and (which come from 

we obtain a bound on the portion from x to X2: 




Xi 


7+1 or^y/S |ari-a;2F 


noting that |x — Xi| = La> 4:Lk, and |xi — X 2 I = Lb> ALk- We claim that 
Lt - cfLt + Lt> {La + 2Lfc + > (|x - X 2 I V ALk)^. 


(2.63) 


(2.64) 


The hrst inequality may be reduced to the minimal case La = Lb = AL^ by comparing 
derivatives of both sides. Then it becomes A"^ — (j) + A"^ > 10^, which is easy to verify 


for -0 = |, 0 = I. The second inequality is clear. Since the inequality is strict in (2.64), 


there is a sufficient leftover power of 7 to control the factors of |x — Xi| and and 
obtain an overall bound 

^(|x-x2|v4L,)’^ . 4"^ (2.65) 

which is sufficient to repeat the argument when summing over (x2, X2) in the next block. 
In this manner we may proceed sequentially through a chain of m blocks. 
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The point z must lie either in one of the blocks or in one of the connecting graphs. 
If a connecting graph is required to pass through 2 ;, then we may take, for example, 
Lb = \xi — z\-\-\z — X 2 \ and the argument goes through, achieving decay along a path that 
includes 2 ;. If 2 ; is in a block, then the path from x to z to y may entail a double traverse 
of the block, but this is covered by the factor of 2 in (2.64). Furthermore, there is also 


the possibility of a self-line, where connects a block to itself. Such graphs contain 
a point z' that i s a distance > L^-i from 5, so \x — z'\ + \z' — y\ > 2Lk-i = Lk, and then 
summing ( |2.60 ) over z' and y gives a bound on the sum of all such graphs. 

Thus we may account for the self-lines with a factor 1.1, say, per block. Another factor 
per block controls the sum over which part of the graph contains z. These considerations 


allow us to obtain from (2.64) the desired decay distance (|a; — y\^-\y — z\) V L^, along 
with a factor 3“^"^ for a chain of m blocks from x to y. Summing this over m > 1 


leads to the desired bound, (2.60) 


follows immediately. This completes the proof. 


The bound (2.61) on the sum of graphs on all scales 

□ 


Corollary 2.8. Under the same assumptions as Theorem\2.'}[ 


IF. 


(k) 


<SbFP(B)\\ < (q7)'-C 


( 2 . 66 ) 


Proof. Graphs contributing to the difference go from x to y via a point z such that 
|x — > Lfc_i, \y — z\> Lk-i. We may bound the norm by estimating the maximum 

absolute row sum of the matrix. This means hxing x and taking the sum over z and y 


of (2.61). A combinatoric factor controls the sums, where L = |x —+ \y — z\ > 

2Lfc_i = Lfc, and (2.66) follows. 


□ 


We also need to control the difference F^\b) — F'^\B) in norm, so that when 
isolated blocks are dehned via the condition dist(specF^^ ^^(F),F) > it is still safe 
to build the random walk expansion for the Schur complement with respect to A. 

Theorem 2.9. With tjj = (f = ^, and e = sufficiently small, assume that |A —F| < 
£fc/2. Then 

||Ff (F) - fP{B)\\ < cJ\X - E\. (2.67) 


yfc)/ 


Proof We have f|^^ = In addition to the explicit ap¬ 
pearance of A, the matrices B^^\ depend on A for k >2. We already have 

control of the graphs contributing to these expressions by Theorem HI An inductive 
argument will allow us to control also the difference when we change A to F. 


We begin by proving an analog of Theorem 2.7 to control the sum of differences of 

' - ' '^(k) 

graphs, i.e. each graph is evaluated at A and at E and the difference taken. Let Sx^y 
denote the sum of the absolute values of all difference multigraphs that contribute to 

We claim that 


S^xl < —■ 2-^|A - F|, 

Sk 


( 2 . 68 ) 
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and hence that 




E\. 


(2.69) 


i<k 


Consider the case k = 1. Redoing the proof of (2.10) for differences, we obtain a sum 
of graphs wherein a difference 

KH^'K - A)-‘ - {WW - = (A - - A)-‘(W'A'> - (2.70) 

appears in place of the corresponding matrix element of — A)“^ or — E)~^. 

Compared to what we had before, there is a combinatoric factor of m for the choice of 
i, a factor |A — E\^ and an extra 1/e from the additional — E)~^. We obtain in 

place of ( 2.62[ ) 


- E\ < - E\, 


x,y — 


\ e J 


2e 


(2.71) 


which verihes (2.68). For the second inequality of (2.71), let L = \x — y\ \/ 1. One can 
check that |(L V 2) > — |, and then the claim may be verihed by counting powers 

of 7. 


For step k, we apply the difference operation to each term of (2.55), and also to 
R(^),C(^) in the expression — X)~^C^^\ Each matrix is 


covered by (2.69), by induction, and this leads to an incremental factor of £ ^|A — E\^ 


comp ared t o before. When we difference the explicit factors of A in (2.55), we obtain 
as in (2.70) a new factor of — £')“^|A — E\. This leads to an incremental factor 

e'/^\X — E\, compared to before, coming from the bound —i7)“^|| < e'/^. Thus in 


all cases, we get no worse than a factor ^|A — E\. This completes the proof of (2.68), 


(2.69). 


(k) (k) 

Note that (2.69) provides an estimate on the matrix elements ol E/ — E/, , so that 


E|[T’--F7i. <cy|A-£|, 


7(^1 


'^1 


(2.72) 


and hence 


IF. 


(k) 


7(^1 


< Cd-\X 

e 


The same bound applies to ||f|^^(F) — E^ 


F|, (2.73) 

*'*'^(5)11, since in this case we are just looking 


at a subset of the collection of multigraphs (the ones that remain within a distance 
< Tfc-i of B). □ 

In order to prove decay of eigenfunctions^ we need a bound similar to Theorem 113 
on the sum of multigraph contributing to , the eigenfunction-generating kern e 
denote the sum of the absolute values of all multigraphs for 
for its dehnition. 


see (2.57),(2.58) 


Le 
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Theorem 2.10. With tp 
|A — i?| < ek/‘2. Then 


= |,0 = and £ = 7'^ sujficiently small, assume that 


g{k) < (1 




■ 2 " 




j<k 


(2.74) 

(2.75) 


Note that along with the graphical terms c ontrol led by this theorem, contains 
an identity matrix at the block for see (2.57). The theorem only controls the 

nontrivial graphs in G 


(fc) 


Proof. We have already established in Theorem 2T the neede d bou nds on G^'^f 
and also which appear in the random walk expansion (2.55) for — A)“^. 

When estimating we have a picture similar to Figureexcept without 


the line. In the proof of Theorem 2.7, gaps in decay and factors ^ at blocks were 


controlled by (2.64). There we had decay on either side of each block, which allowed 
us to prove uniform fractional exponential decay. Here we need to work with just the 
line on one side. Consider the situation in Figure Allowing for combinatoric factors 


[(H^(fc) _ a)-1],5. 



Figure 3: A graph from the eigenfunction-generating kernel. 


as in (2.63), we obtain a bound 

^ (2.76) 

with \x — y\ = La > We claim that 

Lt-^L*>HLa + LtrK (2.77) 

The worst case for this inequality is when La = 4Lfc, where it reduces to 4'^ — 0 > |5^, 
which can easily be verihed for p = ^,'0 = |. Hence we may bound (2.76) by 


(2.78) 


Repeating this argument along a chain of blocks, and then along the chain of operators 
— A) that go into the definition of G'^^\ we obtain ( |2.74[ ). The bound 

(2.75) on the sum over scales follows immediately. □ 
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3 Proofs of Main Results 


We have established probabilistic estimates on resonant blocks, giving decay both in the 
spatial extent and in the volume of the block. Furthermore, we have obtained graphical 
estimates that imply decay of eigenfunctions away from resonant blocks. Let us put 
these tools to work by proving our main theorems. In the hrst section we prove our the¬ 
orem on density of states. In Section 3.2, we introduce an “energy-following procedure” 


that serves to construct every eigenstate through a process of successive approximation. 


In Section 3.3, we prove the needed probabilistic bound on discriminants. In Section 3.4 


we complete the proof of fractional exponential decay of the eigenfunction correlator. 
In Section 3.5[ we extend the method to complete the level-spacing result. 


3.1 Density of States 

Theorem 3.1. There exists a constant b > 0 such that for 7 sufficiently small the 
following is true. For any 0 < 5 < 7^^^, any rectangle A, and any interval I = 
[E — 6/2,E + 5/2], let N{I) denote the number of eigenvalues of H in I. Then 

EN{I) < 4|A|exp(-6|log7|2/^|log5|^/^). (3.1) 


The focus of this theorem is on the behavior as 5 —)■ 0, in particular we see that 
EA^(/) vanishes faster than any power of 1/| log(5|. If one considers instead a “macro¬ 
scopic” interval with 5 > e: = 7'^ = 7^/"^, then the following bound holds: 

EN{I) <4V6\A\. (3.2) 


This estimate may be obtained by a small modihc ation of step 1. Dehne using 
the criterion that tx € [E — 6, E + 6]. Then as in (2.2) we have a bound of 2\/6 on 
the probability that a position is in Weyl’s inequality limits the movement of 

eigenvalues to 0(7), hence N{I) is bounded b y 2|i ?^^^|. The exp ectat ion of is the 
sum on X of the probability that x G R^^\ so (3.2) follows from ( |2.2[). 

Proof of Theorem 3.1 We follow the argument just given for (3.2), but work with 


R^^'^ instead of R^^\ Here k is defined by the condition < 5 < Sk (recall that 

The probability that x G is 


Sk = By assumption, 5 < 7^/^ = e = e\. 

bounded by 


< e^^k-2 < £ 


X^X — 


^2kx _ ,|(2'=’/')1/3 


(3.3) 


using (2.48), q = q/8^, x = V'/S- We may write this as 


exp(-g| log^ll log^^fc+il^/^) < exp(-g| loge:| | log^ 

= exp(—g| logep/^l log5|^'^^). 


(3.4) 


and after summing over x, we obtain (3.1). □ 

Remark 1. The bounds (3.1),(3.2) imply the following unified estimate, valid for 
0 < 5 < 1 : 

EN{I) < 4|A| exp(—6| log(7 V (5)|^'^^| log5|^/^). (3.5) 
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Remark 2. The fractional power of | log5| inside the exponential is a result of the 
discrepancy between x (the probability exponent) and '0 (the graphical exponent). The 
relation y = -0/3 (or something similar) is necessitated ultimately by the ^ exponent 


in (2.18), which leads to the volume limitation in Dehnition 2.1 This leads to limits 


on X from (2.19) and (2.44), depending on how large n is. 


3.2 Energy-Following Procedure 

We will need a modihed procedure that chooses energy windows to contain successive 
approximations to eigenvalues of H. Observe that the spectrum of H is necessarily 
contained in the union of 2(i7-neighborhoods of each G A. (Note that 


< 2(i7, see (1.2),(1.3).) Let us choose one particular x and put Ei equal to one 
of the two eigenvalues of h^- Obviously Ei depends on x and on the choice of 
eigenvalue, but we suppress the dependence in the notation. We may then dehne the 


2.1 


Schur complement 4? , as described in Section 

Recall that blocks were formed by taking connected clusters of sites such that 
dist(ii^i, {ty, —ty}) < e (based on nearest-neighbor connections). Note that if we wish 
to estimate the probability that x, y are contained in the same block, as before there 
must be a sequence of distinct nearest-neighbor positions in connecting x to y. The 

probability that any of these positions (other than x) is in jg bounded by 3 a/£, as 

before. By construction, x is in so there is no small probability associated with 
the event {x G Thus we need to dehne vill, probability that y belongs to 

the block based at x (with a particular choice of Ei taken as hxed). We have that 


< (12dv^)l"-^l. 


(3.6) 


Compared with (2.3), we have one fewer factor of 12dy/e, for the reason just described. 


Graphical bounds such as (2.10) are unaffected by the new procedure. 


To go on with the second step, we will need the truncated version of the Schur 
complement matrix. For |A — i?i| < e/3 we dehne as before F^\b) by restricting the 
sum of graphs that dehnes E^\ including only those that remain within a distance 
< Lq= k oi B. Thus in this step, we consider no graphs of This 


^0 2 
means that 


^^\b) does not actually depend on A. But note that interaction terms 

(b/ 


in V that hop between positions of B are included in F/ (B). Since ||R|| < 2d'y, the 
eigenvalues of E^^\b) are within 2d'f of the “bare” eigenvalues Ety, y ^ B. In a manner 
similar to what was done in the hrst step, we put E 2 equal to one of the eigenvalues of 
Fj^^{B) in [El — e/3, Ei -|- e/3]. There can be no more than n{B) choices for E 2 given 
our initial choices of {x,Ei). We will discuss the counting of choices Ei,E 2 ,... below; 


for now let us observe simply that in view of (2.20), there is a rapid decay of probability 


with n{B), so the choice of E 2 is under control. (There is a potential for redundancy 
here, because an eigenvalue may be reached via multiple initial choices of {x,Ei)] this 
is not a problem as long as every eigenvalue is covered at least once.) 

The process continues in the step in a similar fashion. At this point, there is 
a sequence of choices x,Ei,, E^-i, and the associated increasing sequence of blocks 
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containing x, which may be denoted B^^i,..., We look at the spectrum of 

■^{k—l) I 

A 


Fr ^ (-Ba;,fc-i) in [Ek -1 — Ek-i/S, Ek -1 + £fc_i/3], for A in this same interval. We look 


for solutions of the equation A G specF^^ as this closely approximates the 

(k—1) 

exact equation A G spec E^ ^ that determines eigenvalues of if - see the fundamental 
lemma (we discuss the existence of such solutions below). Theorem 2.9 shows that 


the matrix ^\Bx,k-i) depends only weakly on A: it obeys a Lipschitz condition 


with constant je. By Weyl’s inequality, the same is true of its eigenvalues. Hence 
one may sweep A through the interval above, and all the solutions to the equation 
A G specF^^ ^’{Bx^k-i) in [Ek-i —e:fc_i/3, Ek-i + e:fc_i/3] may be used as the next set of 
approximate eigenvalues. We put Ek equal to one such solution. (Note that by requiring 
\Ek — Ek-i\ < e:fc_i/3, we have that \Ek — Ej\ < e:j/2 for j < k, because the sum of shifts 
Ei/S for j < i < k is less than £j/2. This ensures that we never leave the “safe” zone 
covered by Theorems 2.7, 2.9, and [2.10 ) Each choice may th en se rve as the central 


energy for the next Schur complement E^'^ as defined in Section 2.4[ It is important to 


note that when we shift to Ek, we shift E^'^ —)■ for the random walk expansions 


U) 


at level j < k as well. We have the flexibility to do this because in step j, Theorem 2.7 
controls graphical expansions constructed at any A such that \ \ — Ej\ < Ej/^. 

There are some new aspects to the probability estimates for this “energy-following” 
procedure, as compared to the previous fixed-energy procedure. We need to adjust for 
the fact that by picking Ek close to one of the solutions to A G s]iec F^~^\Bx^k-i) , we 
lose the ability to obtain smallness of the probability of B^ k-ii as is resonant 

to Ek by construction. More precisely, there is no new smallness in the probability 
estimate; when B^^k-i was formed in the previous step, any nontrivial structure was 


associated with smallness of probability (as we saw already in (3.6) in step 1). Note, 


however, that once B^^k-i (the block containing x) is fixed, resonance probabilities 
for the other blocks of involve random variables that are independent of those 

involved in defining Ei,... ,Ek and B^k-i- This is why we always use the truncated 
matrices F^\Bx,j) to define the next set of central energies. Due to this independence, 
the probability estimates for the blocks not containing x (in particular Proposition |2.3[ ) 
are no different from before. Thus the Theorem 2.4 bound (2.22) holds for blocks not 
containing x. 


We prove a modified version of Theorem |2.4| that applies to the block containing x, 
for fixed choices of Ei ,..., Ek. Let B^^k denote the block of x in step k, with B^fi = {x}. 


We introduce an altered version of the inductive definition (2.20) that is suited to the 
energy-following procedure: 


P“>(B) = 


f(3v^) 

l(3v^) 




n(B) 


if X G F; 

if X ^ F 




f)ik)(Tj\ — ) if is isolated on scale k — 2 and x ^ B] 

~ otherwise. 


(3.7) 


As before, n{B) is the number of positions in B, and Bi,... ,Bm' are the subcompo¬ 
nents of B on scale k — 1. It should be clear that this definition eliminates any factor 
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corresponding to an event where x, or a block containing x, is resonant on any scale. As 
explained above, those events hold by construction in the energy-following procedure, 
so no smallness in probability is available. We also need the corresponding weighted 
sum of probability bounds, c.f. ( 2 . 21 [): 




E 




(3,8) 


B containing x and y 
Bx,k\Bx k — 


Theorem 3.2. Given V’ = 3 , IG x = Then for e sufficiently small, and qi = 

qk ^ we have 


Q'tl < c' 


9fc(h-?/|VLfc_2)^/2 


(3.9) 


Note that no smallness of probability can be expected for a block that consists of the 
single site x] smallness arises only when other blocks resonate with an energy from the 
block at X. But we see that for nontrivial blocks, the probability estimate is similar to 
Theorem 2.4, but with a factor of | in the exponents in (3.8),(3.9). Another difference 
is the requirement in (3.8) that \ 7 ^ 0- This is needed because we may 


obtain the minimum decay length Lk _2 only if new resonant blocks have been added to 
B^^k-i in step k. 

Proof of Theorem 3^ We follow the main steps of the proof of Theorem |2.4[ For the 
case k = 1, note that B^^i \Bxfl 7 ^ 0 implies that n > 1. Hence we may counterbalance 
the weighting by probability that n — 1 positions other than x are 

resonant, i.e. for n > 1. Recalling that gi = |, we hnd that the 

sum of all trees with n sites is bounded by Then, since n > |x — 1 /| -|- 1, we 

obtain (3.9) for k = 1. The basis of this argument is that the missing factor of S^/e can 


be absorbed into the bound simply by halving the exponents, as was done in ( |3.8 ),(3.9) 
When we move to the inductive step, we consider a tree graph on blocks, as in Figure 
and repeat the argument beginning with (2.33). Only one of the scale k — 1 blocks 
Bi,... ,Bm' can contain x. We use the induction hypothesis (3.9) if a block contains 


X, and the stronger bound (2.22) from Theorem 2.4 otherwise. However, there is no 


guarantee that Bx,k-i \ Bx,k -2 7 ^ 0- So we insert a partition of unity according to the 


from (3.9): 


< ^qk\x-y\>^/2 ^ 


value of j, the first scale at which B^j = B^^k-i- Then B^j \ B^j-i 7 ^ 0, and we have 

(3.10) 

j<k j<k 

which is valid for x ^ y. (The sum over j can be controlled via the decrease qj —)■ qk in 
the exponent.) Compared with (2.41), the exponent is halved and there is no minimum 
decay length. Note that Bx,k \ B^^k-i 
a “halving argument’’ 


7 ^ 0 implies that m' > 1 , so we may make 


as above to improve the dehcient bound (3.10) on the block 
containing x, at the expense of the others. Specihcally, we use the inequality 


l\xi 


■ I/ll + ^2 - I/2I V Lk-3 > 


- yi\ V Lk-s) + |(|X 2 - 1 / 2 I V Lfc-s). (3.11) 
























Then it shou ld be clear that we can apply the remaining arguments proving bounds 
such as (2.41) on the terms m>2. We obtain an estimate 

m^th Similar arguments will lead to the minimum decay length Lk -2 if m = 1 and 

m' > 1. The condition \ Bx,k-i 7^ 0 eliminates the case m = m' = 1. We obtain 
(3.9) after summin g ove r m. □ 

we may translate this theorem into an estimate on the 


As in Corollary 


2.6 


probability that y belongs to the block of x on scale k. Summing (3.9) over j as in 
(3.10), we obtain the following result: 


Corollary 3.3. Given i/j 
q > 0 such that 


X = '0/3, let e be sufficiently small. Then there is a 


^iV — 


(3.12) 


3.3 Discriminants 


The following result is important for controlling the energy-following procedure, and it 
will also enter into the proof of our main level-spacing theorem for H. It is an analog of 


(2.19) that applies to discriminants, instead of determinants. It will be used to control 


the probability that the level spacing of a block is less than 6. Let us put 


r(u)=disc(Ff*(B) 


(3.13) 


where disc(M) = ~ discriminant of a matrix M with eigenvalues 

Ai,..., Xn- It is well-known that disc(M) is a homogeneous polynomial of degree N{N— 
1 ) in the entries Mij. 

Proposition 3.4. Let B be a block of volume n. Consider the discriminant of the 
Schur complement matrix F^\b). For any 5 > 0, 


P (|r(u)| < ■ A{Anp+\ 


Furthermore, 


P 


(^mjn |Ai -Xj\<S^ < ■ 4(4n)’"+0 

where Ai,..., X 2 n are the eigenvalues of F^\b) . 


(3.14) 


(3.15) 


Proof As explained in the proof of (2.19), the entries of F^’^\b) are equal to those 


of H, plus terms of size /e) that are independent of u = {ux}x&b- Hence r(u) is 

a polynomial of degree < in u. (Here n is the number of positions in B, so that 

N = 2n.) 

We apply the Brudnyi-Ganzburg inequality ( |2.16 ) to r(u). We need to hnd a point in 
parameter space where a lower bound on r(u) can be proven. Take Ui = (1, 2,..., -n,). 
Then the eigenvalues are close to (±a/ 2, ±-\/5, • • •, ±\/l -|- n^). The minimum gap is 
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\/5 — \/2 Ri .82, and allowing for eigenvalne movement (9(7^/e:), we obtain a crnde 
lower bonnd 

r(ui) > (4/5)^"'. (3.16) 


Then, following the argnments for (2.17)-(2.19), we pnt U = [— k = 4n^, and 

a; = {u : r(u) < (3.17) 

Since snp^/ |r(u)| > (d/S)^”^ and \U\ = (2n)”, we obtain 

l/(4n2) 


kn < 


( 5 / 4 )^"'5^34^2 


■4n(2n)’" < ) . 2(2n)^^\ 


(3.18) 


A factor 2"’ converts this into a probability, and we obtain (3.14). Now, restricting 
to the snpport [-1,1] of the random variables Wj, the spectrnm is contained in [—|, |], 
so all eigenvalne diff erence s are bo unded by 3. Hence, if minj<j |Aj — Aj| < 6, then 

implies (3.15) and the proof is complete. □ 

to our energy-following procedure, in the context of a par- 


r(u) < 5234^2 _ (|3.14 


Let us relate this resul 
ticular choice of Ei,..., E^. 

Definition 3.5. A block B of volume n is said to be autoresonant on scale k if 

n < , diam{B) < L^, and 


min |Aj — AjI < ek = , 

i<j 

where Aj, i = 1 ,..., 2n are the eigenvalues of F^\b). 


(3.19) 


Corollary 3.6. Given x, Ei,... ^E^ in the energy-following procedure, the probability 

r tp/z jcy 

that any block containing x is autoresonant on scale k is bounded by for j 

sufficiently small. 


Proof. Given the restrictions on diameter and volume, the number of possible blocks 
B is bounded by 

^2 20 ) 


(2Lfc)'^” < 


For a given B, (3.15) shows that the probability is bounded by 

. 4(4„)"+i < .4(4„)“+i, 


(3.21) 


A combinatoric factor 2"' allows us to £x n. The corollary follows by taking the product 
of these bounds. □ 


3.4 Eigenfunction Correlators 

We wish to prove fractional exponential decay of \g)a{y)Pa{.z)\, which is a strong 

form of localization. One should be able to prove that \g)a{y)Pa{.z) \ is exponentially 
small, except on a set whose probability tends to zero rapidly with \y — z\. One would 
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introduce connectedness conditions on a rapidly growing sequence of length scales as 
in na, for example. Here, we focus on decay of the averaged eigenfunction correlator, 
which in our procedure can decay no more rapidly than do the probability estimates, 
i.e. as a fractional exponential. 

Note that the eigenvalues of H are simple, with probability 1. This can be seen by 
noting that the discriminant is nonzero, for a particular set of well-spaced u’s (of order 
|A|), because the eigenvalues remain separated after turning on parts of H that connect 
different positions. The discriminant is a polynomial in u, so it cannot vanish on a set 
of positive measure without being identically zero. 

Theorem 3.7. Let x = |. There exists a constant q > 0 such that for 7 sufficiently 
small, the eigenfunction correlator satisfies 

(3.22) 


for any rectangle A. 


Proof There are three parts of the proof. In the hrst part, we relate the sum on a 


in (3.22) to the energy-following procedure. We will establish that every eigenvalue Aq 


of H can be reached via the energy-following procedure for at least one set of choices 
for X, El, E2, .... Specihcally, there is a set of choices such that E^ = Aq for all k 
sufficiently large (depending on Aq). Initially, we will allow ourselves the freedom to 
pick from amongst all the sites or blocks that have spectrum close to Aq at each scale. 
Then, we will show that there is at least one good choice of x such that the procedure 
converges to Aq when working with the block containing x at each scale. 

We proceed step by step. In the hrst step, consider the set of positions x such that 
is within e/3 of Aq. Clearly, this set must be nonempty, because otherwise there 
would be no spectrum of Hq within e/3 of Aq, and H = Hq + V with ||I/|| < 2d'y. For 
any of these choices of x, we can choose Ei = -Et^ with \Ei — Ao| < e/3. We then 
construct and E^'^ for |A — Ei\ < e/2. By the fundamental lemma, Aq G specF|^\ 

Now consider what happens after having made choices for xi, Ei, E 2 ,..., Ek such 
that Ej G specand \Ej — Ao| < ej/3 for j < k. Note that \Ej — Ao| < 

ej/3 implies that Aq G specE^/J, by repeated applications of the fundamental lemma. 

Let us truncate —>• The difference has norm < < el, by 

Corollary |2.8[ Therefore, at least one block has spectrum close to Aq, in the sense 
that dist(Ao, spec .F_(^^(i?)) < el. For any such block, we wish to choose Ek+i to be a 
solution in [Aq — ek+i/3, Ao-l-£fc+i/3] to the equation A G spec. f{^^(H). Let us order the 
eigenvalues of as Ai(A) < ■■■ < A„(b)(A). For some p, |Ap(Ao) —Ao| < el. By the 


2.9) and Weyl’s inequality, the interval 


Lipschitz continuity of E^\b) in A (Theorem 
[Ao — 2e:fc_|_i/3, Ao -l- 2ek+i/3] maps contractively into itself under the map A —)■ Ap(A). 
By the contraction mapping principle, there is a solution Ek+i to the equation Ek+i = 
Xp{Ek+i) that satishes \Ek+i — Ao| < 2el < ek+i/3. This may be used to construct 
^(fc+i) , and the procedure continues. 
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For Lfc > diam (A), there can be at most one block, and F^^\b) = Also, 

as explained above, all eigenvalnes are simple, so eventnally there will be exactly one 
solntion to A G specF^^^ in [Aq — £fc+i/3, Aq + £fc+i/3], namely Aq itself. Thns we have 
established convergence of the procednre - at least if one is allowed the freedom to 
pick any of the sites/blocks that have spectrnm close to Aq at each scale. However, we 
would like a somewhat stronger statement: that we can reach Aq starting at some x and 
continuing through a sequence of blocks containing x. To see this, start at scale 
k such that F^\b) = F^'^ and Aq G specF|^\ As we proceed downward in scale, we 
can always conclude that at least one subblock is resonant with Aq (in the sense that 
Ejj^i G specF^.^^^(i?) and |F^j+i — Ao| < The argument is the same as the one 

given above, except that (rather than F^^J) is the matrix that is truncated to 

its subblocks. Then at least one subblock has Ej G specF^. with \Ej — Ej+i\ < Sj. 
Proceeding down to the level of individual positions, we hnd one or more choices of 
X that can be used as a starting point for the energy-following procedure, as claimed 
above. This completes the hrst part of the proof. 

Let us dehne to be the number of eigenvalues of H that can be reached via the 
energy-following procedure starting at x, with a resonant region that includes y. In the 
second part of the proof, we obtain bounds on the expectation of - see (3.29),(3.30) 
below. For each fc, the procedure involves a choice of eigenvalue Ek, which we will sum 
over so as to account for every eigenvalue that can be reached starting at x. Once 
is selected, there is a sum over B^^ki the block containing x (along with an associated 
probability - the block sum is a partition of unity in the probability space). We may 
dehne a “stopping scale” k as the hrst k such that all of the following conditions hold: 


(3.29),(3.30 


1. B^k = Bxki where k is the smallest integer with Lj. > diam(A). In other words, 
Bx^k has reached its maximum extent. 

2 . n{Bx,k) < Lt^". 

3. dieim{Bx,k) < Lk- 

4. Bx,k is not autoresonant on scale k. 


After scale k is reached, there is no further enlargeme nt of Bx,k- Furthermore, conditions 
2-4 imply that the eigenvalues of F^^{Bx^k) obey (3.19), i.e. they have a minimum 
spacing Sk- This means that there is only one choice for Ek+i, Ek+ 2 , ■ ■ ■ (because for 
k > k the intervals [Ek — efc+i/3, Ek + Sk+i/d] can contain no more than one solution 
to the equation Ek+i G F^^^^{Bx^k)-) As a result, we may bound the entirety of 
the choices of Ei, E 2 ,... by [2n{Bx^k)]-- (In step fc, there are no more than 2n{Bx,k) 
eigenvalues of (Bx^k), hence no more than 2n{Bx^k) < 2n(Bx^k) solutions to Ek+i G 

specFg^^^^ {Bx,k)-) For a given set of choices Ei, E 2 ,..., the probability that Bx^k has 
size n and contains y is controlled by Theorem 3.2 Likewise, Proposition 3.4 controls 
the probability of autoresonance. 


32 







As in the proof of Corollary 
that Bx,k = B, i.e. that B is t! 


2.6 


we know that P^^\B) is a bound on the probability 


le block containing x at step k (for a given sequence of 
energies Ei,E 2 ,... for the energy-following procedure). Note that if k is the hrst scale 
where conditions 1-4 all hold, then at least one condition fails when k = k = k — 1. 
Thus we consider four cases, and bound 

EAf.,, < E/VW + EAf<5 + E/V® + EAf^, (3.23) 

where each term represents the case where the corresponding condition breaks down 
for k = k. For case 1, this means that B^^k \ Bx,k-i 7^ 0- Let us assume x ^ y for the 
moment. Then we may estimate 






k=l B containing x and y 


< 


E E 

k=l B containing x and y 

Bx,i^Bx,k—i^0 


(fc- 


k=l k=l 


-'?(h-y|vLfc_2)^/2 


(3.24) 


since for 7 small c.f. (3.8). The last inequality is from Theorem 3.2 (recall 

that Qk \i q as k ^ 00 ). Noting that the exponential growth of L ^_2 with k dominates 
k log k from Stirling’s formula, we have that 


^iV 


(3.25) 


For case 2, with k = k — 1 we bound (2n(P))^ by as above. Then since 

n{B^ we may write 


EivS<E E 


[2n{B)f-^P^^\B) 


k=\ B containing x and y 




(3.26) 


k=l 


j<k 


where we have exploited the decay in n that is built into (3.8) and applied (3.10). For 
case 3, we have that diam(P^ ^) > L^, so B^ ^ must contain a point z with \z — x\ > 
(take z = y a \x — y\ > Lp). Then (after absorbing the factor e^"' into Q as above) 

eN2< E «E E «S+ E «E«« 


k-.L^>\x-y\ j<k 2;|2-x|>L^ 


k-L^<\x-y\ j<k 


< ^9(h-!/|vn^)^/2.5 ^ ^q\x-y\x/3^ 

k=l 


(3.27) 
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Here we have used (3.10) again, and since in either term there is a minimum decay 
length we may use some of that decay to control the sums over 2; and k. Finally, in 
case 4, is autoresonant on scale k, which implies that \x — y\ < and n < 


(here n = Therefore, (2n)^ < exp (A:log(2n)) < exp(A:^). We apply Corollary 

(3.28) 


3.6 to bound the probability, and then 


EiVS< 




fc: Lf,>\x-y\ 

Combining the bounds on the four cases, we obtain 


(3.29) 


\i X = y, the above bounds apply when k > 1. But we also need to consider the case 
fc = 0, which means that the block containing x contains no other points. Then there 
are only two eigenvalues that can be reached, and so 


E W < 3. 


(3.30) 


The third part completes the proof of Theorem |3.7[ We follow the proof of ( |3.29[ ),( |33o| ), 
excep t instead of simply counting eigenvalues, we weight them by \ipa{.y)H^a{.z)\. From 
(2.58) we know that any eigenfunction (p of FT can be written as where i 


IS 


(k) 

an eigenvector of ^ with eigenvalue A. Here we may take k = k\/ k, which ensures 
that Ek is resolved and that F^^{B) = F^^^, so that Ek is an eigenvalue of H. 

Let denote the number of eigenvalues of H that can be reached via the energy¬ 

following procedure starting at x, and whose resonant region includes y, z. Then clearly 


IEiV3.,y,2 ^ E and EiVa,,y,2 < E W,2- 


Taking the geometric mean, we obtain 




(3.31) 


(3.32) 


where we have applied the bounds (3.29),(3.30). 
Theorem 2.10 implies that 

5 ;= E I (G£V'‘')fe)(GtV‘‘>)(^) 


< 








5^ 


|z-ziF 

+ (-' 


W 


3 ;,j/i, 2 i 5 


(3.33) 


(Recall that contains an identity matrix in the block This leads to the 

Kronecker 5 terms in (3.33), as we need to consider separately the cases y G 
y ^ Bx,k, and 2; G Bx^k, z ^ B^^k-) Taking the expectation and inserting (3.32), we 
obtain decay from y to z through the intermediate points yi, x, zi. Thus 


E \ipa{.y)y^a{,z)\ < "^1' 


(3.34) 
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for some q' < q/8 (recall that x = t/'/S- the probability decay is slower than the 
graphical decay). Replacing e with 7 *^, we obtain the corresponding bonnd in (3.22). 
This completes the proof of Theorem |3.7[ □ 

We may prove a rehned version of Theorem |3.7| in which only spectrnm in an interval 
[E — 5/2, E + 5/2] is considered. The resnlting bonnd combines the density of states 
estimate from Theorem |3.1| with the decay estimate jnst completed. 


Corollary 3.8. Let x = |. There exists a constant q > 0 such that for 7 sufficiently 
small, the bound 


o: \a£^[E-S/2,E+5/2\ 


(3.35) 


holds for any 5 < e/2 and any rectangle A. Here k is defined through the inequality 


^fc+i < 25 < 


Proof. Adopting a hybrid approach, we take E to be the cent ral energy in the 


Schnr complement np throngh step k. As in the proo f of Theorem 


the probability that x is in is < c.f (3.3). Then we initiate the energy- 


3.1 


we hnd that 


following procedure, picking x G and the corresponding block ^ of R^^\ We hnd 


a solution to A G spec jf) in [E — ej^/3,E + and continue as before. Now 


when the four cases are considered, c.f. (3.23), there is a minimum decay length T^_ 2 , 
and a corresponding bound 


E A^W < for i = 1,... ,4. 


(3.36) 


If we consider (3.24 


, we have k>k,as all blocks are resonant on scale k. In ( 3.26 ), 


we have k > k a n d k < j < k, which also guarantees a decay of the form (3.36). 


Likewise, in (3.27),(3.28), we have k > k, and a similar estimate holds. This completes 
the proof. □ 


3.5 Level Spacing 

Let us now state our main result on minimum level spacing for H. Let {T^a}Q=i,..., 2 |A| 
denote the eigenvalues of H. 

Theorem 3.9. There exists a constant 6 > 0 such that for 7 sufficiently small, 

p(min\Ea — Egl < d') < |A| exp (—6| log 7 |^/'^| log5|^/'^) , (3.37) 

\ / 

for any rectangle |A| and any 0 < 5 < 7 ^/^/ 4 . 

This is a theorem about microscopic level statistics: we are interested in the behavior 
as 5 —)■ 0. We see that the probability goes to zero faster than any power of 1/| log5|. 
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Proof. We modify the argument in (3.24)-(3.30) in the previous proof. Let us dehne 


Nx{6) to be the number of eigenvalues X oi H that can be reached via the energy¬ 
following procedure starting at x, and which have another eigenvalue within 6 of A. 
Then noting that 

(3.38) 


we see that a bound of the form 


P[mm\E^- E^l < Sj < 


E < exp (—6| logyl^/'^l log51^/“^) 


(3.39) 


implies the theorem. 

Let us dehne k by the inequality 


£k+i <4S < £k. 


(3.40) 


For simplicity we are assuming that S < e/4 so that (3.40) holds for some k > 1. This 


is only a new restriction for moderately small |A|, in any event, since the bound (3.37) 

^-6-^|log|A||^|/|log7|3) 


is nontrivial only for <5 < exp 

As in the proof of Theorem 3.7[ the energy-following procedure generates a sum 
of terms - c.f. (3.24),(3.26),(3.27),(3.28) - that can be used to provide a bound on 


the expected number of eigenvalues A satisfying some condition. Before, we had the 
condition that the resonant region includes x, y. Here, we have the condition that 
another eigenvalue is within 5^/4 of A. Let us break up the sum of terms into part A 
(for which k>k) and part B (for which k<k). (Recall that k is the “stopping scale,” 
i.e. the hrst scale at which conditions 1-4 all hold.) 

Consider hrst the part A sums. Our previous bounds on terms with k > k are 
already small enough for our current purposes. Thus for the part A sums, we may 
ignore the condition that dehnes Nx{£k/4), i.e. that another eigenvalue is within £k/4 


of the ei genva l ue th a 


is r ea ched in the energy-following procedure. Taking y = x m. 


each of (3.24),(3.26[),(3.27), (3.28), and restricting to k> k oi k = k — 1 > k, observe 


that in each case a bound is proven that indicates rapid convergence in k. Noting that 
X = '0/3, the bound is or better in each case (times a prefactor k\ or exp(A;^)). 

We may absorb the prefactor with a decrease in the coefficient in the exponent, and 
then we obtain a bound for the sum of all four cases. 

Consider now part B terms, with k<k. For such terms, all of conditions 1-4 hold 
at step k. In particular, B^^k has reached its maximum extent, and F^^{Bx,k) has 

As explained in the proof 
that can follow E^ in the 


a minimum eigenvalue spacing > - see Dehnition |3.5 

of Theorem |3.7[ there is a unique sequence E^+i 


.E, 




energy-following procedure. Let Aq be the limiting eigenvalue. Then as explained above, 
\Ek — XqI < We have the condition that another eigenvalue A' of H is within 6^/4 of 
Aq. If that is the case, then we have Aq € spec Fxg\ X' G specEy'^ with |Ao — A'| <£k/4. 
We claim that there is some A' with \X' — Ek\ < £k/^ and some level k block B' such 


that A' G spec F^^{B'). This follows from making shifts Fy’ —>■ F^^’ —>■ and 


(k) 


:^P) 


mi 
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applying 


\F? - Fe, 


(fc)| 


< q-|A' - Ek\ < Cd- (efc/4 + 4) , 
e e 


7 


\E 


(k) 


©b-Fe’( s)ll < c: 


(3.41) 


which hold by Theorem 2.9 and Corollary 2.8 Now B' cannot be B^^ki because 
F^^{Bx^k) has minimum level spacing Sk, and Ek G spec {B^^k)■ 

The probability that both B^^k and B' are blocks of is bounded by {Bx,k)P^^"’ {B') 
(By construction, the probabilities and P in (2.20),(3.7) correspond to 

independent events that must occur if B^^k and B' are both blocks of So we may 

sum (B') over all possibilities for B' for each ^ that arises in the energy-following 
pro cedu re. This is bounded by the sum over y of the probability that y G \ B^^k- 


By (3.3), this is bounded by |A|e'^'^fe- 2 . The net result is that the sum of part B terms 


for E,Nx{ek/4:) is bounded by the same four sums (3.24jy^26),(3.27),(3.28), times an 


overall factor of |A|£'^^*- 2 . Taking x = y, we have by (3.30[) that the combined sums 


are bounded by 3, and this leads to an overall bound of 3 |A|£'^^fc -2 on the sum of part 
B terms. 

Adding the above bounds for parts A and B, we obtain 


EN^{S) < + 3\A\e'^^t-2_ 

Let us make a provisional assumption that |A| < P. Then we have 

< + 3e<'"P 

2 -2P/3 _ 1 > 1 guniming over a; G A, we obtain that 


(3.42) 


using 


F(mm |P« - Efd\ < Sk!^ < 


(3.43) 


(3.44) 


The p rovisi onal assumption |A| < £ P is valid in our proof of (3.44), since if it is 


false, (3.44) is automatically true. Recalling that Sk = , we may write 


= g-lg/Sillogelllog^efcA/^ _ g-(q/5)|logep/4|log£fe|i/4^ 


Then since 5 > £fc+i/4 = ef j^-, we have a bound | log(5| < 2| loge^l for e = small. 


(3.45) 


Thus after inserting (3.45) into (3.44), we obtain the statement of Theorem 3.9, □ 


Remark 3. The second term on the right-hand side of (3.42) is subdominant for small 


5, as shown above. If we were to keep it explicit, it would lead to a term proportional 
to |Ap, corresponding to cases where two separate blocks are resonant to each other. If 
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one wishes to obtain a level-spacing bound for 5 < 1, that term becomes the 

dominant one. To see how this affects the level-spacing estimate, look directly at the 
level-spacing condition for {tx}x&K-i and apply (2.2): 

P{mm\tx-ty\ <25) < i|A|(|A| - 1) ■ 2^^ < (3.46) 

x^y 

The eigenvalues of H differ from {P.tx}x&A by 0 ( 7 ), by Weyl’s inequality. Hence 


F( min \Ea — < 5) < V^|A|^, for 7^^^/4 < 5 < 1. 

a^/3 


(3.47) 


Remark 4- We see a decrease —)■ in the exponent of the probability bounds 

when switching from density of states estimates to level spacing estimates. This arises 
because the discriminant is a polynomial of degree 4n^, whereas the determinant has 
degree n. With some optimization of our procedure, one should be able to obtain 
exponents close to respectively, with -0 close to 1. These correspond to the 

worst case in a tradeoff between e” and or respectively. 
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